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INTRODUCTION 

It is a pleasure to participate in celebrating Franco Iachello's 50th birthday. I 
wish him many happy returns. 

At least two physical supersymmetries are known to exist in nature. One, discov- 
ered by Franco and his coworkers, makes connections between different nuclei [1] . The 
other, atomic supersymmetry, is the subject of part of this paper. It can be viewed 
as a symmetry-based approach to the construction of an effective central-potential 
model describing the behavior of the valence electron in atoms and ions. 

This paper begins with some background information and a summary of results 
in atomic supersymmetry. The connection between the supersymmetric Coulomb 
and oscillator problems in arbitrary dimensions is outlined. Next, I treat the issue 
of finding a description of supersymmetry-based quantum-defect theory in terms of 
oscillators. A model with an anharmonic term that yields analytical eigenfunctions is 
introduced to solve this problem in arbitrary dimensions. Finally, I show that geonium 
atoms (particles contained in a Penning trap) offer a realization of a multidimensional 
harmonic oscillator in an idealized limit. The anharmonic theory presented here 
provides a means of modeling the realistic case. 



SUPERSYMMETRIC QUANTUM MECHANICS 

This section provides some background in supersymmetric quantum mechanics 
[2] and establishes notation. 

A quantum-mechanical hamiltonian Hg is said to be supersymmetric if it com- 
mutes with N supersymmetry operators Qj and if it is generated by anticommutators 
according to 

{Qi,Qk} = S jk H s . (l) 

The operators Hs and Qj form the generators of a superalgebra denoted by sqm(AT). 



For the purposes of this paper it suffices to consider the special case N = 2, with 
superalgebra sqm(2). Define the linear combinations 



Q = \[l{Qi + iQ2) , Q ] = yj\ (Qi - 1Q2) ■ (2) 
Then, the supersymmetric hamiltonian can be written 

H s = {Q,Qh . (3) 

For one-dimensional quantum systems, the superalgebra sqm(2) admits a two- 
dimensional representation. Write 

o-Cl 0) • — ft ,°_) • w 

with 

A = - it/72 . (5) 

Here, t/ 7 denotes dU/dx for some function [/ = U(x). In this representation, the 
supersymmetric hamiltonian H$ contains two components h+ and h_, referred to as 
the bosonic and fermionic hamiltonians, respectively. These satisfy the equations 

d 2 

h±q ±n = [-— + V±(x)]q ±n = e n ^± n , (6) 

with 

V ± (x) = (\U'f T \U" . (7) 

Using the above relations, some general properties of the supersymmetric system 
can be obtained. First, for unbroken supersymmetry the ground-state energy is zero. 
Second, except for the ground state (which appears in the spectrum of h+) the bosonic 
and fermionic spectra are degenerate. Finally, the supersymmetry generators Q, Q> 
map degenerate states from the two sectors into each other. 



ATOMIC SUPERSYMMETRY: EXACT LIMIT 



Consider the Schrodinger equation for the hydrogen atom. In spherical polar 
coordinates, the equation separates into angular and radial parts. The angular part 
gives the spherical harmonics, while the radial part can be expressed as 



r d 2 1 1(1 + 1) 1-1 , , 



Here, atomic units are used, and 
y = 2r , E n = 



2n 2 



Xni&r) = rR n i(r) 



The radial wave functions R n i (r) are given by 



Rnl(r) = -2 



rr 



r(ra-Z) 
r(n + / + l)J 
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i) exp (-D L -S(^ 
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(8) 



(9) 



(10) 



In this expression, L^*\x) are the Sonine-Laguerre polynomials, rather than the 
more restricted Laguerre polynomials (for which a must be integer). This distinction 
is important in subsequent sections. The Sonine-Laguerre polynomials are defined by 



v^/ T(n + a + l) 



Ln {X) ^ X) p\(n-p)\T( P +a + l) ■ 

With I fixed, the terms in brackets in Eq. (8) can be reinterpreted in terms of the 
hamiltonian h+ of Eq. (6) as h+ — e n . The requirement that the ground-state energy 
be zero permits the separation of h+ and e n , so that the supersymmetric partner 
hamiltonian h- and the supersymmetry generator Q can be found [3]. These can be 
given by specifying the function U of Eq. (5), which is 

U(y) = j^-2(l + l)\ny . (12) 

Explicitly, the hamiltonian h- looks like h+ but with the constant /(/ + 1) replaced 
with (/ + l)(/ + 2), i.e., 

h _- h+ = w±±l . (13) 

This shift implies that for each / the eigenfunctions of h- are R n j+i, where n > 2. 
Together with the continuum states, they form a complete and ortho normal set. 

This formalism can be given a useful physical interpretation as follows. Con- 
sider the case / = 0. The bosonic sector then describes the s orbitals of hydrogen. 
Since the spectrum of h- is degenerate with that of h+ except for the ground state, 
and since the supersymmetry generator Q acts on the radial part of the hydrogen 
wavefunctions but leaves the spherical harmonics untouched, h- describes a physical 
system that appears hydrogenic but that has the Is orbitals inaccessible. One way of 
realizing this in practice is to fill the Is orbitals with electrons, thereby excluding the 
valence electron by the Pauli principle. The element with filled Is orbitals and one 
valence electron is lithium. This suggests that h- should be interpreted as an effec- 
tive one-body hamiltonian describing the valence electron of lithium when it occupies 
the s orbitals. At this level, the description cannot be exact because most of the 
electron-electron interactions are disregarded. However, these can be introduced as 
supersymmetry-breaking terms. One procedure for this is outlined in a later section 
of this paper. Even in the absence of such terms, some experimental support for this 
atomic supersymmetry can be adduced; see ref. [3]. 

By redefining the energy of the 2s orbital in lithium to be zero, the hamiltonian 
h- becomes a suitable choice for a bosonic hamiltonian of a second supersymmetric 
quantum mechanics. The fermionic partner can be constructed, and an analogous 
interpretation to the one above can be made. This suggests the s orbitals of lithium 
and sodium should also be viewed as supersymmetric partners. The process can be 
repeated for s orbitals and can also be applied for other values of the angular quantum 
number /, leading to supersymmetric connections among atoms and ions across the 
periodic table. In the exact-symmetry limit, these connections all involve integer 
shifts in / and are linked to the Pauli principle. See ref. [3] for more details. 



OSCILLATOR REFORMULATION: EXACT CASE 



Before describing a method for incorporating supersymmetry-breaking effects, it 
is appropriate to discuss an alternative formulation of the exact- symmetry case using 
harmonic oscillators. This section outlines the connections that exist between the 
radial equations for atomic supersymmetry generalized to arbitrary dimensions and 
those for the supersymmetric harmonic oscillator [4]. 

Consider first the <i-dimensional Coulomb problem. Upon separation into an 
angular and a radial part, the radial equation appears: 



d 2 1 , (Z + 7 )(Z + 7 + l) 1 



- + 



dy 2 y r 
Here, atomic units are used, and 



Vdnl(y)=0 . 



(14) 



y = 2r 



E, 



dn 



l=\{d-Z) . 



2(n + 7 ) 2 

The radial wave functions are given by 

v d nl(y) = c dnl y l+ ^ +1 exp(-y/2(n + 7 ))4-S +1) (*/(" + 7)) , 



(15) 



(16) 



where ca n i is a normalization constant. 

As before, this one-variable equation can play the role of the bosonic hamilton- 
aian in a supersymmetric quantum mechanics. Appropriately redefining the energy 
zero so that the ground state has vanishing eigenvalue permits the identification of 
and hence of U, Q and For example, 



U(y) = 



y 



Z + 7 + 1 



2(/ + 7 + l)lny 



(17) 



Just as for the d = 3 case, the fermionic and bosonic hamiltonians differ only by the 
replacement of / by / + 1, so that 



h. 



h+ = 



2(/ + 7 + l) 



(18) 



As required, all the results of the previous section are recovered if d = 3, i.e., 7 = 0. 

Consider next the supersymmetric harmonic oscillator. For convenience, anal- 
ogous variables to the Coulomb-problem quantities y, d, n, /, h± are now denoted 
by upper case symbols Y, D, N, L, H±. Thus, upon separation of the angular and 
radial parts of the Schrodinger equation for the D-dimensional harmonic oscillator, 
the radial equation is obtained as: 



d 2 2 (L + r)(L + r + i) 

dY 2 Y 2 



2E DN jv DNL (Y) = . 



(19) 



Here, atomic units have again been used for simplicity, and the oscillator is assumed 
to have unit frequency. The radial variable is now Y, and 



E DN = l(2N + 2T + 3) , r=(D-3)/2 



(20) 



The radial wave functions are given by 



V DNL (Y) = C DNL Y L + r +i exp(-y 2 /2)L^t 1 / 2 2) (^ 2 ) , (21) 

with Cdnl a normalization constant. Note that the usual expressions for the har- 
monic oscillator in three dimensions are recovered when T = 0. 

If Eq. (19) is used to define the hamiltonian H + of a supersymmetric quantum 
mechanics (a redefinition of the energy zero is again needed), then the supersymmetry 
is specified by a function U given by 

u(Y) = y 2 -2(L + r + i)i n y . (22) 

It then follows that H- differs from H + by the replacement of L with L+l. Therefore, 

H _ _ H+ = !(£±£±11 . m 

So far, four eigenspectra associated with the supersymmetric Coulomb and oscil- 
lator problems have been introduced, defined by the four hamiltonians h + , h-, H + , 
and H_. The hamiltonians h + and h- are related by the map / — > / + 1, and the 
hamiltonians H + and H_ are related by L — > L + 1. The next step is to relate the 
(i-dimensional Coulomb problem to the D-dimensional oscillator, i.e., connect h+ to 
H+. 

It can be shown [4] that an eigenfunction of h + can be transformed by a one- 
parameter mapping into an eigenfunction of H + . Explicitly, the functions Vdni and 
Vdnl are connected by the equation 

v dn i{(n + ~f)Y 2 ) =K DNL Y 1 ' 2 V DNL {Y) , (24) 

where Kdnl is a proportionality constant and 

D = 2d-2-2X , N = 2n-2 + \ , L = 21 + A . (25) 

The integer A is the mapping parameter. Notice in particular that only oscillators in 
even dimensions appear. Note also that if it is desired that more than one eigenfunc- 
tion of h + be mapped into the eigenspace of a specified H + , then the possible choices 
of D, N, L, and A can become constrained. 

The existence of the one-parameter map between h + and H + combined with the 
supersymmetry maps evidently establishes connections between any two of the four 
hamiltonians h +: H +: and H_. More details can be found in ref. [4]. 



BROKEN SUPERSYMMETRY AND QUANTUM-DEFECT THEORY 

As noted above, atomic supersymmetry in the exact limit is not physically real- 
ized because the valence electron interacts with the core electrons by more than the 
Pauli principle. This section discusses the incorporation of supersymmetry-breaking 
effects in the context of alkali-metal atoms. 

One important effect of the interactions between the valence electron and the core 
is the change in energy eigenvalues relative to the hydrogenic case. In alkali-metal 



atoms, the Rydberg series [5] provides a simple formula for the measured energies, 
given by 

*«- = -2^ ■ (26) 

In this expression, 

n* = n-6(n,l) , (27) 

where 5(n, I) is called the quantum defect. For a fixed value of / and increasing n, it 
turns out that the quantum defects rapidly attain asymptotic values: 5(nJ) ~ 5(1). 

The changes in the energy eigenvalues imply that the exact atomic supersym- 
metry is broken. The breaking can be viewed as an additional contribution Hb to 
the supersymmetric hamiltonian Hs of Eq. (4). For example, if h+ arises from the 
radial equation for hydrogen and h- is interpreted as the radial equation for the va- 
lence electron of lithium in the exact-super symmetry limit, then the hamiltonian H 
describing the two systems including supersymmetry-breaking effects can be taken as 

H = H S + H B , (28) 

where Hb has the form 

M° w) (29) 

and Vb(u) is such as to generate the observed energy eigenspectrum of lithium. 

The determination of a suitable Vb is not straightforward. However, it turns out 
that a functional form for Vb can be found that yields analytical eigenfunctions as 
solutions to the Schrodinger equation [6]. It is 

v ( v **(** + !)-*(* + !) , n 2 -n* 2 

Here, /* is a modified angular quantum number given by 

I* = I + - 6(1) , (31) 

where i(l) is an integer parameter shifting the angular quantum number in a manner 
characteristic of supersymmetry. (If desired, 5(1) could be replaced by 6(n, I).) This 
model effectively replaces the hydrogenic radial equation with one of similar form but 
involving n* and /* rather than n and /. 

By construction, the energy eigenvalues are those of the physical atom. The 
resulting eigenfunctions _R*, z *(r) are analytical and are given by 

F(n*-l*) ]h(2r\i\ r\_ (2I . +1) /2r\ 



K*l*(r)- n * 2 [ r(n * + z * + 1) J { n *) ex P(, n *) L n-i-i-i{ n *) ■ ( 32 ) 



The Sonine-Laguerre polynomials enter again because 

n *_/*_l =n -i-i(i)-i (33) 

remains integer. For asymptotic quantum defects 5(1) and including the continuum 
states, these eigenfunctions form an orthonormal and complete set. 

More details about this construction can be found in ref. [6]. (A connection to 
parastatistics is elucidated in ref. [7].) There is a reasonable body of evidence to 



support the notion that the analytical eigenfunctions provide a good model for the 
valence electron, especially in alkali-metal atoms. For instance, transition probabil- 
ities calculated with the analytical eigenfunctions agree with experiment and with 
accepted values [8]. (Some recursion formulae for matrix elements are given in ref. 
[9].) Transition probabilities for other elements, notably alkaline-earth ions, have 
also been obtained in this way [10]. Moreover, these analytical eigenfunctions have 
been used as trial wavefunctions in detailed atomic calculations [11]. Stark maps 
for the alkali-metal atoms can also be calculated using the model [12]. The resulting 
clear anticrossings and small-field quadratic Stark effects for the s and p orbitals are in 
agreement with experiment. For example, the model yields Stark maps for the n = 15 
lines of lithium and sodium that are indistinguishable from the numerical and experi- 
mental results of ref. [13]. (Ref. [12] also studied other possible quantum-mechanical 
supersymmetries involving hydrogen. In particular, a double sqm(2) appears when 
the separation is carried out in parabolic coordinates.) The model is expected to 
break down at short distances from the nucleus (of order of the core size), but despite 
this some dominant features of the fine structure in alkali-metal atoms are correctly 
reproduced and the Lande semiempirical formula naturally appears [14]. 



OSCILLATOR REFORMULATION: BROKEN CASE 

This section presents a reformulation of the analytical quantum-defect model 
in terms of oscillators with radial equation modified by an anharmonic term. For 
generality, the connection between the two is treated in arbitrary dimensions. The 
oscillator models that appear have analytical solutions. The link between the two 
theories is via a three-parameter map. In the limit of vanishing quantum defect, this 
map provides a generalization of the ones of ref. [4] discussed above. For example, 
it can be used to connect a Coulomb problem to an anharmonic oscillator with odd 
dimensionality. 

The first step is to construct the ^-dimensional extension of the quantum-defect 
model of ref. [6]. This is done by adding to the hamiltonian (14) an extra term V^(y) 
generalizing Vb in Eq. (30), with 

v d ( . _ (/*+7)(/*+7 + l)-(/ + 7)(/ + 7 + l) , (n + 7 ) 2 -(n*+ 7 ) 2 m , 
V bW- ~2 + 4(n + 7 ) 2 (n* + 7 ) 2 ' ( ' 

Here, n* and /* are modified quantum numbers given by 

n*=n-5(dj) , I* = l + i(d,l)-8(d,l) , (35) 
in analogy with Eqs. (27) and (31), with 

7=s(d-3) (36) 
as before. The ensuing radial equation in atomic units is 



#_ _ i + (t + tKt + t + i) _ i_ Ednt 

dy 2 y y 2 2 n 



J dn*l 



(y) = o , (37) 



where ^ 

^ = "2(n*+ 7 )2 • (38) 



The radial wave functions solving Eq. (37) are given by 

v* dn *l*(y) = 4^* 2 / r+ ^ 1 exp(-y/2(n* + 7 ))4 2 i;^:+ 1) ( ?/ /(n* + 7 )) , (39) 

where c^ n * z * is a normalization constant. 

To identify a mapping between the quantum-defect theory and an oscillator-type 
model, a (supersymmetry-breaking) term Vg(Y) modifying the harmonic-oscillator 
radial hamiltonian (19) is needed. A suitable choice is 

vjPtn = < £ ' + r »< L ' + r + v- < L + r >< £ + r + x > + 2 ( N - m . (40) 

In this expression, 

r = £(£>- 3) (41) 
as before, and the modified quantum numbers A* and L* are given by 

A* = A — 2A(D, A, L) , L* = L + 2I(D, L) - 2A(D, A, L) , (42) 

where the integer 2/(1), L) is a supersymmetry-type shift and 2A represents a quan- 
tum anharmonicity (which can be viewed as an oscillator 'defect'). The factors of 
two are introduced for notational simplicity in what follows. The extra term (40) 
introduces an anharmonic piece into the oscillator potential, which in turn changes 
the energy eigenspectrum. In atomic units with a unit- frequency oscillator, the an- 
harmonic radial equation becomes 

'-■^ + Y 2 + + T)i y 2 +T + 1) - 2E DN *] VZ N . L .{Y) = , (43) 
where the energy eigenvalues are shifted according to 

E DN * = |(2A*+2r + 3) . (44) 
The eigensolutions for this anharmonic oscillator are 

vW(y) = c* DN . L . r^+ r+1 ex P (-yV2)< 2 + _^ , (45) 

where Cp N * L * is a normalization constant. 

The map between the radial equations for the (/-dimensional quantum-defect the- 
ory and the .D-dimensional anharmonic oscillator connects the eigensolutions v^ n *i* 
and Vp N * L „. It is given by 

v* dnH *{{n* + 1 )Y 2 )=K* DN * L *Y 1 / 2 VZ )N * L *{Y) , (46) 

where Kp N * L ,. is a proportionality constant and 

D = 2d-2-2A, A = 2n+2(A-5)-2+A , L = 2l+2(A—6)—2(I—i)+\ .(47) 

For fixed d, n, /, 6, and i there are three quantities that effectively act as mapping 
parameters: A, A, and /. For the eigenfunctions (45) to exist / must be integer, so 
the supersymmetry-type shift 21 in L is an even integer. Since D must also be an 
integer, A is integer or half-integer. The quantum numbers A and L are also integer, 



which implies that 2(A — 5) + A must be integer. Note that this generalizes the 
exact-symmetry case: the requirement that A be a whole integer is no longer needed 
because half-integral values can be absorbed in the difference 2(A — 8). This means, 
for example, that when the Coulomb problem is treated in the exact limit (S = 0) it 
is now possible to map it into a (modified) oscillator in an odd number of dimensions, 
provided A is quarter-integer valued. As in the exact case, for maps between sets of 
states with specified d and D, further restrictions on the possible values of D, N, L, 
/, A and A may appear. 



SUPERSYMMETRIC OSCILLATORS AND THE PENNING TRAP 

In this section, I demonstrate that geonium atoms provide a physical realization 
of a D > 1 super symmetric harmonic oscillator. For simplicity, the specific case 
D = 2 is considered, although under suitable conditions higher values of D may ap- 
pear. In practical situations the supersymmetry is broken for reasons to be described. 
The analytical anharmonic oscillator model introduced in the previous section should 
provide a good approximation to the exact wavefunctions for this case. Space limi- 
tations prevent more than a sketch of the relevant physics being given here; details 
will appear elsewhere. 

Geonium atoms are formed by a set of charged particles bound in a Penning 
trap [15], which is a suitable combination of a homogeneous magnetic field and an 
electrostatic quadrupole potential. The simplest geonium atom has just one trapped 
particle of charge e and mass m [16]. Successively adding further electrons in the 
trap generates elements of the geonium periodic table. 

For simplicity, consider the idealized Penning trap with electromagnetic fields 
specified in cylindrical coordinates (p, 9, z) by 

B = Bz , ^I-j^-Ip 2 ) . (48) 

The quantity d is a measure of the trap dimension and is to be specified in terms 
of the configuration of the quadrupole electrodes. For (stable) trapping, eV > 0. 
The quantum-mechanical motion of the particle in the field B is that of a harmonic 
oscillator with binding frequency equal to the cyclotron frequency, given in SI units 
by 

\eB\ 

u c = ^ . (49) 
m 

(In fact, there are two oscillators involved in this motion, but only one enters the 
quantum hamiltonian.) Similarly, the electrostatic field generates an axial harmonic 
motion independent of the cyclotron motion, with axial frequency 



eV 

(50) 



d 2 



These are the motions of primary interest here. 

The simultaneous presence of electric and magnetic fields also generates another 
(unbound) circular motion, called the magnetron motion, with frequency u m . For 
simplicity, this is largely disregarded here. The eigenvalue spectrum of the system 
is split by all these interactions and also (for particles with spin S) by the spin 
interaction — S • B. The latter splitting implies the existence of a supersymmetry of 



the type discussed in refs. [17]. This supersymmetry is not directly relevant to the 
discussion here, and the spin degree of freedom is neglected in what follows. 

The combination of the cyclotron and axial motions forms a system of two one- 
dimensional oscillators. This becomes a physical realization of a D = 2 harmonic 
oscillator when the applied electromagnetic fields are chosen such that u c = u z , i.e., 

mV = \e\B 2 d 2 . (51) 

The quantum problem can then be separated in polar coordinates and the radial 
equation has the general form of Eq. (19) with T = — |, and with suitable constant 
factors inserted to allow for non-unit binding frequency and for SI units. 

For fixed L, this oscillator can be used as the bosonic partner H+ in a super sym- 
metric quantum mechanics. The partner hamiltonian H_ is specified by Eq. (23). It 
represents a system having an eigenspectrum degenerate with the bosonic sector but 
with the ground state missing. As in atomic supersymmetry, one practical realiza- 
tion of this is to fill the ground state with particles and invoke the Pauli principle. If 
L = 0, for example, H + describes the S orbitals of the simplest geonium atom with 
one trapped particle. Then, if_ can be interpreted as an effective theory describing 
the behavior of the 'valence' particle in the S orbitals of a more complex geonium 
atom in which the IS (and 2P) orbitals are filled. This interpretation invokes the 
approximation in which particle interactions other than those implied by the Pauli 
principle are disregarded. All the atomic supersymmetries of ref. [3] have analogues 
in this system. For example, there are connections between pairs of geonium atoms 
throughout the geonium periodic table. 

There exists a mapping between the d = 3 Coulomb problem and the D = 2 
harmonic oscillator, as discussed above. The connection between the eigenfunctions 
given in Eq. (24) therefore establishes a correspondence between eigenfunctions of 
elements in the usual periodic table and the geonium periodic table. The map (24) 
is fixed here by setting A = 1, so that 

N = 2n - 1 , L = 21 + 1 . (52) 

Moreover, in the exact supersymmetry limit this map induces other maps involving 
the supersymmetric partners so that all four hamiltonians are linked. 

The supersymmetries are broken by the interaction of the valence particle with 
the 'core' of the geonium atom. These interactions will shift the eigenenergy of the 
valence particle from En = N + 1 (in the level with quantum numbers N, L) to some 
other energy 

E N * = N* + 1 , (53) 

where by definition 

N* = N - 2A(A, L) . (54) 

A model incorporating these exact new eigenenergies and yielding analytical solutions 
has been introduced in the previous section. In the present case, it is obtained by 
adding an anharmonic term to the oscillator hamiltonian, giving the radial equation 
(43) with r = — | and with suitable dimension-correcting factors inserted. The 
analytical solutions are given by a corresponding modification of Eq. (45). 

It is physically plausible to conjecture that the quantum anharmonicity rapidly 
approaches an asymptotic value as N becomes large, i.e., 



A(N,L) ~ A(L) 



(55) 



In this case, the eigensolutions form a complete and orthogonal set. The previous 
section also provides a mapping between this theory and the analytical quantum- 
defect theory for ordinary atoms. Combined with the magnetron and spin splittings, 
the anharmonic model is likely to provide a simple method for calculations of physical 
properties of the valence particle in geonium atoms. 
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